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ABSTRACT 

We construct multimonopole solutions containing N — 1 distinct fundamental monopoles 
in SU (N) gauge theory. When the gauge symmetry is spontaneously broken to [/(l)^ 1 , 
the monopoles are all massive, and we show that the fields can be written in terms of 
elementary function for all values of the monopole positions and phases. In the limit 
of unbroken U{\) x SU{N — 2) x U{1) symmetry, the configuration can be viewed as 
containing a pair of massive monopoles, each carrying both U{1) and SU{N — 2) mag- 
netic charges, together with N — 3 massless monopoles that condense into a cloud of 
non-Abelian fields. We obtain explicit expressions for the fields in the latter case and 
use these to analyze the properties of the non-Abelian cloud. 
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1 Introduction 



The massive magnetic monopoles of spontaneously broken gauge theories have long been the objects 
of considerable study. Although these monopoles arise as spatially extended solutions of the classical 
field equations, they correspond to single-particle states of the full quantum theory. Their dynamics, 
at least at low energies, can be described by a small number of degrees of freedom, just like that of 
the elementary particles of the theory jlj. Indeed, these magnetically charged states can be regarded 
as the counterparts, related by an exact duality symmetry in certain supersymmetric theories, of 
the massive electrically charged states built from the elementary excitations of the theory || . 

Recent studies of low-energy monopole dynamics have shown that when the unbroken gauge 
group has a non-Abelian component there are degrees of freedom that can be attributed to the 
presence of massless non-Abelian monopoles; these can be seen as the dual counterparts of the 
massless elementary gauge bosons ||. In contrast with the massive monopoles, these massless 
monopoles cannot be exhibited as isolated classical solutions, but can be studied classically only 
as part of multimonopole configurations. In the simplest example, an SO(5) solution Q with 
one massive and one massless monopole, the massless monopole is manifested as a spherically 
symmetric "cloud" of non-Abelian fields surrounding the massive monopole. In this paper we 
examine a somewhat more complex class of configurations containing two massive and N — 3 
massless monopoles in the Bogomolny-Prasad-Sommerfield (BPS) @j limit of an SU(N) gauge 
theory. By obtaining explicit analytic expressions for the gauge and Higgs fields we can see clearly 
how the massless monopoles condense into a non-Abelian cloud and can verify that the properties 
of this cloud inferred from the form of the moduli space metric are indeed present. 

The origin of these massless monopoles can be understood by recalling that the magnetically 
charged BPS solutions for an arbitrary gauge group G can all be analyzed in terms of fundamental 
monopoles of various types. The simplest case is when the adjoint Higgs field breaks a group of 
rank r maximally, to U{l) r . There are then r quantized topological charges, one for each of the 
unbroken U(l) factors. Corresponding to each of these is a fundamental monopole solution carrying 
a single unit of topological charge. Each of these is described by four collective coordinates, three 
specifying its position and one corresponding to an overall U(l) phase. These solutions can be 
realized explicitly by embedding the unit SU(2) monopole using a preferred set of simple roots. 
All higher charged solutions may be regarded as multimonopole solutions containing appropriate 
numbers of the various fundamental monopoles. Not only does one find that the energy of the 
solution is the sum of the masses of the component monopoles, but an index theory analysis shows 
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that the number of zero modes, and hence of collective coordinates, is precisely four times the 
number of component monopoles Q. 

This can be illustrated in a particularly simple fashion when the gauge group is SU(N). By 
means of a gauge transformation, the asymptotic value of the Higgs field in any fixed direction can 
always be brought into the form 

$ = diag (t N ,t N -!, (1.1) 

with t\ < ti < . . . < tjv- We may write the asymptotic magnetic field in the same direction as 

Fij = eijkrkQui/r 3 , with 

4-7T 

Qm = — diag (njv^njv-i - n N _ 2 , ■ ■ ■ ,n 2 - m, -m) . (1.2) 

e 

The generalized topological quantization condition implies that the rij must be integers @. 

If the eigenvalues of the asymptotic Higgs field are all distinct, then the unbroken gauge group 
is U(l) r and the rij are the topological charges. Apart from a constant Higgs field contribution, the 
fcth fundamental monopole, with rij = 5jk, is obtained by embedding the SU(2) monopole solution 
(rescaled appropriately so as to give the correct Higgs expectation value) in the 2x2 block at the 
intersections of the (N — k)th and (N + 1 — k)th rows and columns. The resulting SU(N) monopole 
then has a mass 

M k = 27T{tk+1 ~ tk) (1.3) 

e 

where e is the gauge coupling. Although there are other possible SU{2) embeddings, both the mass 
formulas and the zero mode counting indicate that these are merely special cases of multimonopole 
solutions. 

Varying the asymptotic Higgs field so that two or more of its eigenvalues are equal enlarges 
the unbroken symmetry group so that some of the U(l) factors are replaced by a non-Abelian 
group K. The magnetic charge must still be of the form of Eq. Ql,2| ), but the rij that correspond 
to roots of K are no longer topologically conserved charges; in fact, they are not even gauge- 
invariant. According to Eq. (1.3), the corresponding monopoles should become massless in this 



limit. From the point of view of electric-magnetic duality this seems quite reasonable, since one 
would expect the massless gauge bosons carrying electric-type charges in the subgroup K to have 
massless counterparts carrying magnetic charges. On the other hand, one would not expect to find 
zero energy solitons. Indeed, the classical one-monopole solution tends toward the vacuum solution 
as the limit of unbroken symmetry is approached. However, examination || of the moduli space 
Lagrangian that describes the low-energy dynamics of a collection of BPS monopoles suggests that 
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the degrees of freedom corresponding to these monopoles can survive even in the massless limit. 
Specifically, if a number of massless monopoles are combined with one or more massive monopoles 
to give an n-monopole solution whose total magnetic charge is invariant^] under K, the dimension 
of the moduli space, and hence the number of collective coordinates, remains 4n even in the limit 
of non-Abelian unbroken symmetry Q. Furthermore, examination of specific examples suggests 
that the moduli space metric, and hence the Lagrangian, behaves smoothly in this limit. 

As noted above, the simplest examples, containing one massive and one massless monopole, 
arise in the context of an 50(5) gauge theory spontaneously broken to SU(2) x U(l) [|j. These 
contain a massive monopole core surrounded by a spherically symmetric "non-Abelian cloud", of 
arbitrary radius, that can be viewed as the remnant of the massless monopole. Within the cloud 
there is a Coulomb magnetic field corresponding to a magnetic charge with components lying both 
in the unbroken U(l) and in the unbroken SU(2), while outside the cloud only the U(l) Coulomb 
field is present. The solution is described by eight collective coordinates. Four of these are readily 
identified as the position and U(l) phase of the massive monopole. The other four coordinates 
describe the cloud, with three determining its overall SU{2) orientation and one specifying its 
radius. 

An obvious step toward gaining further understanding of these massless monopoles and their 
associated non-Abelian clouds would be to investigate solutions containing larger numbers of 
monopoles. Solutions corresponding to one massless and two identical massive monopoles have 
been studied in SU(3) broken to SU(2) x U(l) fK| |ll| and in Sp(4) = 50(5), also broken to 
SU (2) x U(l) |12}| . (In the latter case the unbroken SU(2) x U(l) is a different subgroup than that 
considered in Ref. In both cases the moduli space metric was found explicitly, but analytic 
expressions for the gauge and Higgs fields could only be found for special configurations. 

The complexity of these solutions is perhaps not surprising if one recalls the rather nontrivial 



form of the Atiyah-Hitchin metric for the moduli space of two identical SU{2) monopoles 13]. 



By contrast, the moduli space metric for two [14], or even an arbitrary number pSl , of distinct 
fundamental monopoles is relatively simple, and so one might expect the corresponding solutions 
for the fields to be more tractable. As we will see, this is indeed the case. We consider here 
SU (N) solutions comprising N — 1 distinct fundamental monopoles. For both the case of maximal 
symmetry breaking, where all of the monopoles are massive, and the case where the unbroken 
group is U (1) x SU(N — 2) x U(l), where all but two of the monopoles become massless, the gauge 

One encounters a number of pathologies when dealing with configurations whose total magnetic charge has a 
non-Abelian component || . Requiring that the total charge be purely Abelian is not really a significant restriction, 
since the additional monopoles needed to meet this requirement can be placed arbitrarily far from the original ones. 
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and Higgs fields can be expressed in terms of elementary functions for all values of the collective 
coordinates. 



We use Nahm's method to construct these solutions [16|. In Sec. II we review the details of this 
construction for BPS monopoles in an SU(N) theory. The implementation of the construction for 
configurations containing many distinct fundamental monopoles is described in Sec. III. In Sec. IV 
we consider the case where the unbroken group is U(l) x SU(N — 2) x £7(1) and obtain explicit 
expressions for the fields. These expressions simplify considerably in the regions outside the cores 
of the massive monopoles. We discuss these asymptotic forms in Sec. V. Section VI contains some 
concluding remarks. There is an Appendix containing details of some of the calculations. 

2 The Nahm construction 

The fundamental elements in Nahm's construction |l(| of the BPS monopole solutions are a triplet 
of matrices T a (t), the Nahm data, that satisfy a set of nonlinear ordinary differential equations. 
These T a then define a linear differential equation for a second set of matrices, v(t, r), from which the 
fields <3?(r) and A(r) can be constructed. In this section we review the details of this construction 
for the case of an SU(N) theory jl7| with the asymptotic Higgs field and magnetic charge given 
by Eqs. ([O]) and (|1.2|) , respectively^. Here, and for the remainder of the paper, we set the gauge 
coupling equal to unity. 

The matrices T a {t) are defined for t± < t < tjy. The tj divide this range into N — 1 intervals. 
On the jth interval, tj < t < ij+i, we define k(t) = rij and require that the T a have dimension 
k(t) x k(t). In addition, whenever two adjacent intervals have the same value for k(t), there are 
three matrices ctj, of dimension k(tj) x k(tj), defined at the interval boundary tj. These matrices 
satisfy the Nahm equation 

= \e ahc [T b , T c ] + £(<*,■).*(* - tj) . (2.1) 
j 

where the sum in the last term (and similar sums in later equations) should be understood to run 
only over those values of j such that rij = rij-x- (The T a are singular at tj+% if |nj+i — rij\ > 2. 
Because we will not be considering such situations here, we will not describe the requirements 
obeyed by these singularities.) 

Having found the Nahm data, the next step is to find a 2k(t) x N matrix function v(t, r) and 
2 We will in general follow the notation of Ref. mM . 



4 



iV-component row vectors Sj(r) obeying the differential equation 



d_ 

"dt 



+ (T a + r a ) <g> a a 



+ Y j a)S J 5{t-t j ) 



and the normalization condition 



1 = J dtv^v + Y, S ] S i- 



Here dj is a 2/c(ij)-component row vector obeying 



OjCij = ctj ■ a — i(aj)oI 



with (oj)o chosen so that the above matrix has rank 1. 
Finally, the spacetime fields are given by 



$ = J dttv ] v + tjSjSj 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



A = -i J dtvWv -i^sfiSj 

j 



dt 



v'X7v — Viv v 



(2.6) 



where the second equality in Eq. Q2.6] ) is obtained with the aid of the normalization condition 
Eq. (p.3[). These satisfy the self-dual BPS equations 



where 



B = D$ 



B a = ^e abc F bc = ^e abc (d b A c - d c A b + i[A b , A c \) 



(2.7) 

(2.8) 
(2.9) 



Equations (|2.2| ) and (|2.3| ) do not completely determine v and the Sj. Given any solution of 
these equations, a second solution can be obtained by multiplication on the right by an N x N 
unitary matrix function of r; this corresponds to an ordinary gauge transformation. In addition, 
there is also some freedom to multiply v and the Sj on the left, with corresponding transformations 
on the Nahm data. Such transformations have no effect on the spacetime fields, but can be used 
to simplify the intermediate calculations, as we will see in Sec. 4. 
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3 Construction of (1,1,..., 1) monopole solutions 

We will be concerned in this paper with solutions consisting of N— 1 distinct fundamental monopoles. 
The rij are then all equal to unity and k(t) = 1 for the entire range of t. Since k(t) is unchanged at 
each of the intermediate tj, there is an a.j and an Sj for each value of j from 2 through N — 1. The 
commutator term vanishes, and so Eq. is easily solved to give the piecewise constant solution 



T(i) = -Xj , tj <t< t j+1 . (3.1) 

The x a have a natural interpretation as the positions of the individual monopoles. The a,j of 
Eq. ( |2.4D are simply two-component row vectors, which we take to be 



2|xj - Xj-i] (cos(0/2)e-^ /2 ,sin(0/2)e^ /2 ) (3.2) 

where 9 and <fi specify the direction of the vector cxj = Xj_i — Xj. 

The next step is to find a2xJV matrix v(t) and a set of iV-component row vectors S k (k = 
2, 3, . . . , N — 1) that satisfy Eq. (|2.2| ). To this end, we first define for each interval tk < t < tt+i a 
function fk(t), with 

f^t) = e (*-*2)(r- Xl ).<T 

h(i) = e^-^-^h^), k>l. (3.3) 

These have been defined so that their values at the endpoints of the intervals satisfy 

fk(t k ) = fk-i(t k ) = g k , k = 2,3,...,N -1. (3.4) 

An arbitrary solution of Eq. ( |2.2D can then be written in the form 

v(t) = fk(t)r] k , h <t < t k+1 , (3.5) 

with discontinuities at the intermediate t k obeying 

r] k = ?? fc -i +g k 1 a [ k S k . (3.6) 

The normalization condition, Eq. (|2.3| ), becomes 

N-\ N-l 

I= H S j S i + E Vk^N km (3.7) 

3=2 k=l 

where 

N k = / dtfl(t)f k (t). (3.8) 



We will find it convenient to distinguish between the first two and the last (N — 2) columns of 
v and the S k , labeling the former by Roman superscripts from the beginning of the alphabet and 
the latter by Greek superscripts that run from 3 to N. We choose the v a to be continuous, so that 

S£ = 0, a = 1,2. (3.9) 

A properly normalized solution for the v a is then obtained by taking 

rjg = JV-VV , a = 1,2, (3.10) 

where 

N-l 

N=Y, N k (3-11) 

fc=i 

and the a are the two-component objects 8 1 = (1,0)* and 9 2 = (0, 1)*. 

Orthogonality of each of the last N — 2 columns of v with the first two implies that 

N-l 

0=^N k r^, M = 3,4,...iV. (3.12) 

k=l 

Together with the discontinuity Eq. ( |3.6| ), this gives 

N-l N-l 

< = ^ E E cmNug^alSf (3.13) 

k=l 1=2 

where 

f 1, 3>l>k, 
c jkl = l-l, j<l<k, (3.14) 
[ 0, otherwise . 

Substituting these solutions for the into the orthogonality condition, Eq. ( |3.7| ), gives 

JV-l N-l 
= E + "'U-K'^j = E ^KijS* (3.15) 

i,j=2 i,j=2 

where 

AT-l 

M 4J = ]T cuiCknjg^NjN-^kN^NngJ 1 . (3.16) 

fe,Z,m=l 

Hence, 

5J = (if-Va)^^ (3.17) 

where f7 is any (N — 2) x (iV — 2) unitary matrixrL The freedom to choose U corresponds to a 
U(N — 2) subgroup of the SU(N) gauge symmetry; the remaining gauge symmetry has already 

been fixed by our choices for the first two columns of v. 

In applying this equation, one must be careful to take into account our convention that the upper index on S 
runs from 3 to iV, while the lower runs from 2 to N — 1. 
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Substituting this expression for the S% into Eq. (gig) gives the rjj and thus, through Eq. (|3.5|) , 
determines v(i). It is then a straightforward, although tedious, matter to substitute these results 
into Eqs. ( |2.5| ) and ( |2.6| ) and thus obtain the fields A(r) and $(r). We will not carry this out 
explicitly for the case of maximal symmetry breaking. However, we note that it is clear from the 
above equations that the result can be expressed in terms of elementary functions. 

4 Solutions for SU(N) U(l) x SU(N - 2) x U(l) 

Our main interest is in the case where the middle N — 2 eigenvalues of the asymptotic Higgs field 
are all equal, so that the unbroken gauge group is U (1) x SU (N — 2) x U(l). If we adjust the Higgs 
field in this fashion, then, as was argued in Ref. ||, the (N — l)-monopole solutions of the previous 
section can be viewed as being composed of two massive and N — 3 massless monopoles, with the 
latter condensing into a non-Abelian "cloud". The massive monopoles are located at and Xjv_ i; 
without any loss of generality we may take these to lie on the z-axis, with zat_i = z\ + R > z\. 
The locations of the massless monopoles are less well-defined. Extrapolating from the maximally 
broken case, one would take these to be the points X2, X3, . . . , xtv-2- However, as we will now show, 
many different choices for these points yield the same solution. 

To begin, note that of the N — 1 intervals into which the range of t was divided, only the 
leftmost (t\ < t < £2) an d rightmost fa = ijv— 1 ^ t ^ tw) now have nonzero width; we shall use 
subscripts L and R to label quantities related to these two intervals. Hence, of the N — 1 integrals 
defined by Eq. Q3.8|), only N% = Nl and Njy^i = Nr are nonzero. Also, the = fk(tk) are all 



equal to unity. As a result, Eq. (3.15) simplifies to 



i,3 

where 



^ = ^ s**t \ 6 . . + a . Ma r\ s » (41) 



M=(N£ 1 + N R 1 ) (4.2) 



is Hermitian. Once a set of Sj satisfying Eq. (fOj) has been found, the and hence v(t), can be 
found from Eq. ( 3.13| ), which now reduces to 



N^M^S*. (4.3) 

3 

With maximally broken symmetry, the monopole positions enter both through the functions 
fk(t) and through the various aj. In the present case, where the middle intervals have zero width, 



the Xfc associated with the massless monopoles enter only through the aj, which, as can be seen 
from Eqs. Q4.1|) and ( f4,3|) , appear only in the combination Ylj a j^j- With this in mind, consider 
two sets of monopole positions and with identical locations for the massive monopoles, but 
with the massless monopoles constrained only by the requirement that dj = Wj^a^, with W some 
(N — 2) x (N — 2) unitary matrix. If Sj gives a solution of Eq. ( [4.1| ) for the former set of positions, 
j = WjkSk is a solution for the transformed set. (Note that this transformation has no effect on 
either v or on the fields themselves.) 

The possibility of performing such transformations implies that the positions of the massless 
monopoles are not all physically meaningful quantities. In fact, these yield only a single physical 

parameter, which can be identified by noting that these transformations leave invariant the quantity 

AT-l 

a ) a i = [°i ' a ~ ia isA = ( x i ~ x »-i) " a + Y l x i ~ x J"-il ■ ( 4 - 4 ) 
j j i=2 

The first term on the right hand side is fixed by the positions of the two massive monopoles. The 
second term is just the sum of the distances between successive monopoles. It is precisely this sum 
that was identified in Ref. |||, from the properties of the moduli space, as the unique gauge-invariant 
quantity characterizing the non-Abelian cloud. It will be convenient to express it in terms of the 
"cloud parameter" b > 0, defined by the equation 

N-l 

2b + R= J2 l x i - x i-il • ( 4 - 5 ) 
i=2 

It is not hard to show that any two sets of massless monopole positions corresponding to the 

same value of b can be transformed into one another. We may therefore define a canonical set of 

positions by placing one of the massless monopoles on the z-axis at Z2 = z\ — b and the remaining 

N — 3 massless monopoles on top of the massive monopole at £jv— 1> The only nonvanishing aj are 

ci2 and (Z3, which can be combined as the rows of a real 2x2 matrix 

. fV2b \ .... 

A = n /o7m tT\ ■ ( 4 - 6 ) 



y/2{b + R). 

With this canonical choice of monopole positions, the solution for the Sj~ given in Eq. ( p. 17 ) 
takes on a particularly simple form. If these N — 2 row vectors are combined to form an (N — 2)xN 
matrix S, the solution corresponding to the choice U = I can be written as 

S=(° n 2 r°.) («) 



. /at_ 4 , 

where the columns (rows) have been grouped in blocks of 2, 2, and N — 4 (2 and N — 4) and the 
2x2 matrix 

S = S^ = (I + AMA)~ 1/2 . (4.8) 
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Equations ( |3.5| ) and ( ff-3 ) then give i>(t), which in a similar block notation takes the form 

(M)(x-»,-Wl)»,o), tl <t<t 2 , 

\f!t-i(t)(N- 1 ' 2 ,N?WL) 1 ",6) , h<t<t N . 
Here N = Nl + Nr and we have defined 

L = MA 2 {I + AM AY 1 . (4.10) 

Note that, as we will see explicitly below, A, M, and L are all diagonal and hence commute. 

These results can now be combined to yield the gauge and Higgs fields. To express these, it is 
convenient to introduce the integrals 

H L = f 2 dt[h{t),Vh{t)} 
Jtx 

K L = l 2 dtth{tf -t 2 N L (4.11) 

Jtx 

and the similarly defined quantities Kr and Hr involving /jv_i(t). The fields can then be written 

■N- 1 ! 2 0\ / aW a® 0\ / N' 1 / 2 0\ 

(LM)V 2 a(3) f a( 2 ) (ML) 1 /* (4.12) 

o o o/ \ o o o/ V o o o/ 

/iv- 1 ^ o O\/0 (1) (3) 0\ /JV-i/2 o 0\ 
$ = t 2 J + (LM)V2 o 0( 3 ) f 0( 2 ) (ML) 1 / 2 (4.13) 

V o o 0/ \ o o 0/ V o 0/ 

where the a*- a ) and 4>^ are the 2x2 matrices 

a« = -^H L + Hi? ) + l[ivV2 )ViV i/2- 

0« = (Kl + ^/0 

a( 2 ) = -l^H^ + JV^H^ 1 ) 

<t>® = N^KlN^ +N R l K R N R l 

a (3) = '-[(VNu-Hr)^ 1 -(VNl-Ul)^ 1 } 

^ = KrN^-KlNZ 1 . (4.14) 

To proceed further we need explicit expressions for the various integrals that we have defined. 
These are most easily expressed in terms of the vectors 

YL = r - xi 

y R = r - xat_i (4.15) 
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and the quantities 



sl = (*2 — *i)]yi| = (*2-*i)yz 
sr = (t N -t 2 )\y R \ = (t N -t 2 )y R . 



Straightforward integration yields 



N L = y L 1 smhs L e~ SL y L '°' 



N R = y R 1 smhs R e SRfR - (T 
K L = ^V L X YL ■ o- 
Kr = -Vr^Yr ■ cr 

-Wl 2 (yL x <r) (sinh s L cosh s L - s L 



(*2 ~ h)e~ 2s ^ (T - N L 
(t N -t 2 )e 2s ^ a -N R 



WR 2 (yR x a) (sinh s R cosh s R - s R ) . 



Substituting these expressions into Eq. ( 4.14| ) yields the remarkably simple formulas 

1 



a (3) 



(V x a) 
-\CV. a) 

'-Fa 
2 

If 

2 



where 



V = y L 



coth sl 



SL 



sinh sl 



+ YR 



COth Sr 



SR 



sinh 



and 



F = YL - - Yr - + (t 2 - ti) (1 - y L ■ cr coth s L ) + (t N - t 2 ) (1 + y R ■ a coth s R ) 

VL VR 

Combining the expressions for Nl and N R with Eq. fl4.2j), and using the fact that 



(YL - Yr) ■ cr = Ra 3 



we obtain 



M = (y L coth SL + yRCothsR + Ra 3 ) . 

Combining this with Eqs. ( |4.6[ ) and ( [4.10| ) yields 

_ 2b + R- Ra 3 

yL coth SL + yR coth s R + .R + 26 ' 

Note that the cloud parameter 6 enters the solution only through the matrix L. 
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5 Asymptotic behavior 



In this section we will examine in some detail the solutions that we have found. From the form of 
Eqs. (|4.12j ) and ( [4.13| ), it is clear that for any iV > 5 all SU (N) solutions are essentially embeddingsQ 
of solutions with one massless and two massive monopoles in SU (4) broken to U(l) x SU (2) x U(l). 
Therefore, without any loss of generality we can simplify our notation by specializing to the case 
N = 4. Each adjoint representation elementary multiplet of the theory can then be decomposed 
into five massless fields (an SU(2) triplet and two singlets) together with a pair of massive doublets 
and a pair of massive singlets. Because of the "twisting" of the topologically nontrivial Higgs field, 
this decomposition will not in general have a simple correspondence with the matrix components of 
the fields. However, one might hope that matters would simplify in the region outside the cores of 
the massive monopoles (i.e., the region where sl and sr are both much greater than unity), where 
the massive fields would be expected to be exponentially small. 

The first hint of this simplification comes from noting that Nl and Nr each have one exponen- 
tially large eigenvalue. Hence, the eigenvalues of N = Nl + Nr are in general both exponentially 
large, implying that those of iV -1 / 2 are both exponentially small. The only exception to this oc- 
curs when the large-eigenvalue eigenvectors of Nl and Nr are almost parallel to each other, which 
happens only near the line joining the centers of the two massive monopoles. If we exclude the 
region close to this intermonopole axis, no elements of the matrices F, M, and L are ever expo- 
nentially large. We then immediately see from Eq. ( 4. 14] ) that the terms containing a^ 3 ) and (f>^ 
are exponentially small. 

If we therefore restrict ourselves to the region where sl,sr 3> 1 (so that we are outside the 
massive cores) and 

"' "" ~ R - ^>e~ 2sL +e~ 2sR (5.1) 



R 

(to avoid the intermonopole axis), we may approximate the fields by the block diagonal form 

A« 
A( 2 ) 



(5.2) 



•-( o *">)• < 5 - 3 > 

With the fields written in this form, their group theoretic interpretation is fairly clear. The traceless 

parts of the nonvanishing blocks correspond to two commuting SU{2) subgroups, one of which (the 

lower right, as we shall see) must be the unbroken SU(2). The massless U(l) fields are then 
4 For the special case R = and ti = — tjv, these actually reduce to embeddings of the SO(5) solution of Ref. W 
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contained in the traceless part of the other block and in the two traces. Examining Eqs. ( |4.12j ) and 
( [4.13| ) and recalling that the dependence on the cloud parameter is only through L, we see that 
A^ 2 ) and $( 2 ) may be fe-dependent, but and are not. (Because the fields are traceless, the 
6-dependence must be entirely in the traceless parts of A^ 2 ) and & 2 \) 

Because of the factors of N~ 1 / 2 , the analysis required to obtain the asymptotic form for the 
upper left block is in general somewhat tedious. However, the calculation simplifies considerably 
if e 2sL /ul is either much less than or much greater than e 2sR /ur (which is the case in almost of 
space.) If we define the unit vector n to be equal to y R (yl) m the former (latter) region, then, 
up to exponentially small corrections, the gauge and Higgs fields are 

A(1) = 2[(yL V +y R T-R 2 ] m X 9r) + X iH + ' <T ^ + lf Vfl ' cr ' fl '^ (5 - 4) 



and 



2yR) \ 2 J \ 2y L J 
while the asymptotic field strength is 



(5.5) 



(In fact, it is not hard to show that anc [ B^ 1 ) must be of this form (although with a more 
complicated expression for n) at all points outside the massive cores.) 

Although the calculation of the fields in the lower right corner is in principle straightforward, 
somewhat lengthy manipulations needed to put the result in a simple form. We leave the details 
of these to the Appendix, and state the results here. From Eqs. ( [4.12 ) and (4.18) we obtain 



A< 2 ) = -i(LM) 1 / 2 V x a (ML) 1 / 2 (5.7) 
where we can now use the asymptotic forms 

V = y L + y R (5.8) 

M = (y L + y R + Ras)- 1 (5.9) 



and 



L = 2 l +R ;!?Z h - H2b + R - Ra 3 ) . (5.10) 
VL + VR + R + 2b 



For the Higgs field we obtain the particularly simple form 



$( 2 ) = t 2 I + L 1/2 



1 /l — q-tr\ 1 / 1 + q • a 



2y R \ 2 J 2y L 



L 1 ' 2 (5.11) 
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where q is a unit vector with components 

(yL)a + (yR)a 

'la ' 



1,2, 



(5.12) 



Vl-Vr „ — o 

r ' a — 6. 

(Note that at large distances q, like n, approaches the radial unit vector.) The field strength is 



B (2) = L 1 / 2 



2yl + y L 



l + qcr 



yr h (jL + yR) 



2y 



VR 



1 — q • a 



h 



-f L 1 ' 2 



^VhVR 

(5.13) 



1,2, 



fa 



(5.14) 



where h is defined by Eq. (|5.10| ) and 

V(VL + Vr) 2 ~ R 1 <?a 

(VL + y/?) 03 - R, a = 3. 
It is instructive to consider several limiting cases. First, suppose that 6 = 0. From Eq. (4.23) we 
see that L is then proportional to 1 — 03 so that, except for a constant term in the third rows and 
columns of all fields vanish. The solutions are then essentially embeddings of SU(3) —*U(l)x U(l) 
solutions I IS]. Since the unbroken subgroup of SU(3) is Abelian, it should be possible to choose a 
gauge in which the asymptotic fields are simply superpositions of single monopole fields. Indeed, 
these can be written in the form 





$ = l/f J (r) 



1 

2y R 













B = C/f 1 (r) 



/ Yr 






Yl 



YR 
2^ 













*2 

tx 
\ 






Ui(r) 



(5.15) 



V 



CAW 



(5.16) 





[Here, in order to make the U(l) x SU(2) x C/(l) structure of the theory clearer, we have reordered 



the rows and columns to correspond with the order in Eqs. ( |l.l|) and (|l.2j ). With this reordering, 
the unbroken SU{2) is contained in the middle 2x2 block, and a fundamental monopole solution 
with a single nonzero rij corresponds to an embedding in a pair of adjacent rows and columns.] 
Viewed as an SU(3) solution, this corresponds to a configuration containing one each of the two 
massive fundamental monopoles of the theory, located at points xi and X3. Viewed as an SU(4) 
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solution, it can be interpreted as containing a massive fundamental monopole with rij = 5ji at xi 
and a superposition of a massive monopole with rij = 5j3 and a massless monopole^] with rij = 5j 2 
at X3. Even though the underlying SU(3) solution has purely Abelian long-range fields, the long- 
range part of the SU(4) solution is non-Abelian in the sense that the unbroken SU{2) subgroup 
acts nontrivially on A^ 2 ) and $( 2 ). The crucial point, however, is that the SU(2) orientations of 
the Coulomb field centered at xi and the Coulomb field centered at X3 are aligned, so that the net 
SU(2) component is a purely dipole field that falls as R/y 3 for y = (y^ + yn)/2 3> R. 

Next, consider the case 6 3> R. Now L is approximately proportional to a unit matrix, with 

2b 



0(R/b) , 



(5.17) 



yL + yn + 2b 

while h ~ 1/(26 + yi + Thus, in the region where y^ and yn are much less than b the fields 



can be written as 



$ = m 












t 2 + 



2y R 



V 












* 2 ~ 2k 



U 2 (r) + 



(5.18) 



B 



1 2^ 








2^5 








\ 



\ 






U 2 (r) + 



(5.19) 



where the dots represent terms that are suppressed by powers of R/b, yr,/b, or yn/b. These are 
just the fields expected for two massive monopoles, with topological charges rij = 5j\ and rij = 8js, 
centered at xi and X3, respectively. In contrast with the previous case, the SU(2) components 
of their two magnetic charges are not aligned, and so the unbroken SU{2) contains two Coulomb 

fields rather than the dipole field of Eq. ( 5,16| ). These non-Abelian Coulomb fields disappear when 

5 Viewing the massless monopole as being centered at X3 is a gauge choice; moving it to xi simply corresponds to 
a re-ordering of the rows and columns. 
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y S> b, where we obtain 



<S> = U 3 \r) 



(u-± o o o \ 

t 2 o 

t 2 o 
V o 



U 3 (r) + 0(b/y 2 



(5.20) 



o t x + l y ) 



/A 

/ 2^ 



B = C/.7 1 fr) 





V 



\ 








U 3 (r) + 0(b/y 3 ). 



(5.21) 



Thus, at distances large compared to b the magnetic fields exhibit the behavior that would be 
expected if a single massless monopole were added to the two massive fundamental monopoles. 
However, this massless monopole is not manifested as a localized structure with a well-defined 
center. Instead, we simply have a transition from a "cloud" region of size ~ b containing non- 
Abelian Coulomb magnetic fields to an outer region where these fields are cancelled. 

In this discussion we have excluded the region close to the intermonopole axis. To explore the 
fields in this region, one must go back to the equations of the previous section. [Equations (|5.8| - 
[5.10 ) are not valid approximations here, even outside the monopole cores.] Doing so, we find that 
the fields do not have the simple block diagonal structure of Eqs. fl5.2| ) and (|5.3| ). In addition, 
some components of A become exponentially large as one approaches the axis, while $ turns out 
to be rapidly varying. However, these are essentially artifacts of our choice of gauge. The net 
magnetic charge of our solutions is a unit charge in the unbroken U(l) that is contained in A^. 
The long-range twisting of the Higgs field must then be topologically equivalent to an embedding 
of the SU(2) hedgehog configuration in the corresponding SU(2); this can be seen in the behavior 
of the unit vector n that appears in Eqs. (p.^|- [5.6[ ). However, this cannot be the whole story. 
Because the two monopoles have different £7(1) charges, there must be some additional twisting of 
the Higgs field near each of the monopoles. The conventions that we have adopted are such that 



this inevitable additional twisting is confined to the narrow region near the axis where Eq. ( 5.1 



does not hold. In order that D<£ and B not become large, these rapid variations in the direction 
of <3? must be compensated by large values of A. That this actually happens can be verified by 
evaluating the field strength along the axis. One finds that, up to exponentially small corrections, 
the magnetic field along this axis is independent of b and can be put in the form of Eq. ( |5.16| ). 
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6 Concluding Remarks 



In this paper we have shown how the Nahm construction can be used to obtain explicit multi- 
monopole solutions corresponding to N — 1 distinct fundamental monopoles in SU(N). In the case 
where the symmetry is broken maximally, to [/(l)^ -1 , these solutions are described by 3(N — 1) 
gauge-invariant parameters that specify the positions of the component monopoles; these combine 
with N — 1 overall U(l) phases to give the full set of collective coordinates. Even though these 
solutions have no rotational symmetry at all, the gauge and Higgs fields can be expressed in terms 
of elementary (i.e., rational and hyperbolic) functions for arbitrary values of these parameters. 
Of particular interest is the behavior of these solutions in the limit where the unbroken group 
is enlarged to U(l) x SU(N — 2) x U(l) and TV — 3 of the component monopoles become mass- 
less. Examining the solutions, one sees no trace of the individual massless monopoles, but only 
a single "non-Abelian cloud". Indeed, most of the massless monopole coordinates become redun- 
dant in this limit, with the solutions depending only on a sum of intermonopole distances that is 
conveniently described by the cloud parameter b. Not only are the massless monopole positions 
somewhat ill-defined, but so, in a sense, is their number. As we have seen, the SU(N) solution 
that nominally contains N — 3 massless monopoles is essentially equivalent to an embedded 5t7(4) 
solution containing a single massless monopole. 

For values of b that are large compared to the distance R between the massive monopoles, 
the cloud is rather similar to that found in the 50(5) case. The non-Abelian Coulomb magnetic 
fields due to the two massive monopoles are present inside the cloud, but are screened at distances 
much larger than b. On the other hand, when b <C R these Coulomb fields are present outside the 
cloud (i.e., at points more than a distance b from either of the massive cores), but are aligned with 
each other in such a manner that the entire solution is simply an embedding of a purely Abelian 
configuration. 

These non-Abelian clouds and their properties clearly warrant further study. For the configura- 
tions we have considered here, those parts of the solution associated with the cloud are particularly 
simple. The dependence on the cloud parameter is contained in a single matrix L, while corre- 
sponding components of the Higgs and gauge fields depend on the same function of the spatial 
variables and differ only in their tensor structure. In addition, the actual functional forms outside 
the massive cores are relatively simple. This simplicity suggests that it might be feasible to de- 
termine the structure of the non-Abelian clouds associated with more complex configurations, at 
least in the regions outside the cores of the massive monopoles. The most compelling questions are 
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associated with the conjectured electric-magnetic duality. From this point of view, the massless 
monopoles are clearly the duals of the massless non-Abelian gauge bosons and their superpartners. 
How do these particles reflect the strange properties of the massless monopoles? 

We are grateful to Choonkyu Lee, Kimyeong Lee and Changhai Lu for helpful conversations. 
We thank the Aspen Center for Physics, where this work was begun. This work was supported in 
part by the U.S. Department of Energy and the U.S. National Science Foundation. 



Appendix 



In this Appendix we outline the manipulations that lead to Eqs. fl5.ll ) and ( 5.13| ) for <I>( 2 ) and 



B^ 2 - 1 . Throughout, V, M, and L should be understood to be given by the asymptotic forms in 
EqB. |Tfjrijj ). 



The first step is to obtain some useful identities. With the aid of the law of cosines, we obtain 

V* = 2(l + y L -y R )= {yL + m)2 - R2 (A.1) 

VhVR 

Because we have chosen the massive monopoles to lie along the z-axis, the third component of any 
vector can be obtained by noting that Rw% = w • (y^ — yr). In particular, many of the subsequent 
results make use of the identities 

\ (vl + vr) (vl + vr) 2 - R 2 i a 9 x 

iVLh = R 2RY L (A - 2) 



and 



, (vl + Vr) {VL + VRf-R 2 / a n\ 

{mh = R— + 2kV R (A - 3) 



Next, we define a vector f by the equation 

M 1/2 aM 1/2 = 



(VL + VR? ~ R 2 f (A.4) 



Explicit expressions for the components of f are given in Eq. ( |5.14[ ). The inner products of this 
equation with y^ and y# then lead, with the aid of Eqs. (|A.2|) and (|A.3|), to 



M l ' 2 y L ■ aM 1 ' 2 = ( S^Z±1) _ M ( A .5) 
V 2y L J 

and 

M^ 2 y R ■ aM'/ 2 = (^^) + M (A.6) 

where q is defined by Eq. ( |5.12| ). Combining these last two equations with the results of Sec. 4 
gives the asymptotic expression for $( 2 ) shown in Eq. (5.11). 
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To obtain an expression for the magnetic field we need the derivatives of L and M. Making use 
of the fact that these are both diagonal matrices, we find that 



L -i/2 VL i/2 = .Ilvl- 1 = -- 

2 2 



hV 



and 



M~ 1/2 VM 1/2 = -^MVM" 1 = -^MV 



where h = yi + y# + R + 26. In addition, we note that 

L = 1 — hivr 1 

Using the last three identities, we can decompose the magnetic field B^ 2 ) as 



(A.7) 



(Ai 



(A.9) 



^abc 



9 fe 4 2) + i4 2) 4 2) 



L l/2 \ M l/%A) M l/2 + hM y%(B) M l/2] L l/2 



(A.10) 



where 



^abc 
1 



-^(V x a) c - ^{M-V^Af 1 / 2 , (V x <r) c } + i(V x <r) 6 M(V x er) c 



2 MV b - a b d b V a ) + \v b {M, (V a a b - V b a a )} + % -{\ x a) b MV a a h 



(A.ll) 



and 



^ahc 



^V b (V x a) c - x a) b (Y x cr) 



F V • £T - V 2 a a 



(A.12) 



The last equation, together with Eqs. (A.4- A6| ), immediately gives the terms in Eq. ( |5.13| ) propor- 
tional to h. To obtain the /i-independent terms, it is useful to write 



One then finds that 



M 



b^ 



VL + VR 



Vl + Vr 



VLVRy 2 {vl + Vr) 2 - R 2 VlVrY 2 



g°3 



Oil) a 



(m)a 



2yL 2vr 2 



(A.13) 



(A.14) 



which leads to 



M l/2 6 (^) M l/2 



{VL)a ( q • cr + 1 \ (y R ) a /q • cr - 1\ 1 



+ 



2y 2 



R 



J+-M 
2 7 2 



VR VL 



a3) 

(A.15) 



After some algebra, one finds that the quantity in square brackets multiplying M vanishes, thus 
leading to the final expression for B( 2 ), Eq. ( gig ). 
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